Excitons in one-dimensional Mott insulators 
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We employ dynamical density-matrix renormalization group (DDMRG) and field-theory methods 
to determine the frequency-dependent optical conductivity in one-dimensional extended, half-filled 
Hubbard models. The field-theory approach is applicable to the regime of "small" Mott gaps which 
is the most difficult to access by DDMRG. For very large Mott gaps the DDMRG recovers analytical 
results obtained previously by means of strong-coupling techniques. We focus on exciton formation 
at energies below the onset of the absorption continuum. As a consequence of spin-charge separation, 
these Mott-Hubbard excitons are bound states of spinless, charged excitations ( "holon-antiholon" 
pairs). We also determine exciton binding energies and sizes. In contrast to simple band insulators, 
we observe that excitons exist in the Mott-insulating phase only for a sufficiently strong intersite 
Coulomb repulsion. Furthermore, our results show that the exciton binding energy and size are not 
related in a simple way to the strength of the Coulomb interaction. 

PACS numbers; 71.10.Fd, 71.35Cc, 72.80.Sk 



I. INTRODUCTION 

Excitons in conventional hand insulators are well- 
described by Wannier theoryld. A Wannier exciton is 
a charge neutral optical excitation made of an elec- 
tron in the conduction band and a hole in the valence 
band, bound together by the Coulomb attraction be- 
tween them. In inorganic semiconductors like GaAs the 
typical binding energy, as defined by the energy difference 
between the exciton and the band edge of the particle- 
hole continuum, is several meV. This should be compared 
to the band gap itself, which is of the order of one eV. 
Concomitantly the typical size of a Wannier exciton is of 
the order of 100 A, which is almost two orders of mag- 
nitude larger than the lattice spacing. We note that al- 
though the total spin of an optically excited exciton is 
necessarily zero, it is composed of two quasiparticles car- 
rying spin- 1/2. 

In quasi one-dimensional materials like conjugated 
polymera3 the situation is quite different. Here the 
electron-electron interaction accounts, to a substantial 
degree, for the optical gap itselfl as well as for the for- 
mation of excitons. The exciton binding energy in, e^™ 
polydiacetylenes is found to be of the order of 0.5 cVO'EI 
and is thus comparable to the optical gap of 2.4 eV in 
SBCMtUpolydiacetylene chains diluted in tbeir monomer 
matrixala. The exciton size was estimatedtl to be 12 A 
and is thus comparable to the length of the unit cell 
of 5 A. These facts suggest that electron-electron inter- 
actions will play an important role in any theoretical de- 
scription of excitons in these materials. 

Realistic models for conjugated polymers must account 
for the effects of both electron-electron and electron- 
phonon interactions. The interplay between these makes 
the reliable calculation of the optical conductivity a very 
demanding task. As a first step it is therefore nat- 
ural to investigate the effects of the two mechanisms 
separately^. The optical conductivity for models with 



only electron-lattice coupling such as the celebrated Su- 
SchrieffeipHeeger model has been widely analyzed in the 
literatureQ and is well established. As an extension of this 
approach, electron-electroninteractions have been taken 
into account perturbativelya. In contrast there are com- 
paratively few reliable results for the optical spectrum in 
models that talce account of the sizable electron-electron 
interaction&~E3. The interaction-drives these systems into 
a Mott-insulating ground stateE3. The scarcity of results 
is due to the difficulties associated with the calculation 
of excited state properties in one-dimensional Mott insu- 
lators. Therefore, it is of considerable interest to develop 
reliable methods for the investigation of optical excita- 
tions in correlated electron systems, and to determine 
the optical conductivity of one-dimensional Mott insula- 
tors. 

In this paper we focus on the calculation of the opti- 
cal conductivity in models with electron-electron interac- 
tions only. We study an extended Hubbard model with 
nearest and next-nearest neighbor density-density inter- 
actions; the model is further specified in Sec. ||. 

We employ two recently developed numerical and ana- 
lytical techniques to determine the real part of the optical 
conductivity over the full frequency range and analyze 
exciton properties without suffering from finite-size limi- 
tations; for a first application to the Hubbard model, see 
Rcf. [ij. In Sec. Ill we first test the dynamical density- 
matrix renormalization group (DDMRG) by applying it 
to the limit p£ large Mott gaps where analytical results 
are availableEj. We obtain excellent agreement between 
numerical and analytical results, and confirm the clear 
and simple physical picture of an exciton as a bound 
state of a double occupancy and an empty lattice site in 
a background of singly occupied sites. We then move on 
to the generic case of intermediate Mott gaps and find 
qualitatively the same physical behavior. 

In the regime of small Mott gaps, finite-size effects and 
finite resolution of the DDMRG start to hamper the nu- 
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merical analysis. Therefore, in Sec. [V, we carry out a 
weak-coupling field-theory analysis of the problem. Us- 
ing the form-factor bootstrap approach we determine the 
optical conductivity in the field-theory limit. Here, the 
spin sector does not couple to the current operator so that 
it is sufficient to analyze the charge sector only. The ex- 
citon is then described as a bound state of a holon and 
an antiholon, which are the elementary charge excita- 
tions in the theory. The resulting picture for small Mott 
gaps remains very similar to the cases of intermediate to 
large gaps. We even find quantitative agreement between 
field theory and DDMRG results for intermediate Mott 
gaps where the applicability of field theory is not a priori 
expected. 

In Sec. we discuss two fundamental properties of 
Mott-Hubbard excitons, their binding energy and size, 
in greater detail. In contrast to Wannier excitons in band 
insulators, Mott-Hubbard excitons exist only when the 
intersite Coulomb repulsion exceeds a certain threshold. 
In general, the exciton binding energy is not related in a 
simple way to the strength of the Coulomb interaction. 
We analyze a new correlation function which allows to 
define the size of an exciton in correlated electron sys- 
tems. Our analysis provides a comprehensive picture of 
excitons in one-dimensional Mott insulators. 



In our conclusions (Sec. VI) we address implications of 
our results for the theory of 7r-conjugated polymers. 



a spin-flip symmetry and a spatial-reflection symmetry 
(through the lattice center). Therefore, each eigenstate 
of (|l|) has a well-defined parity under charge conjugation 
(Pc = ±1) and spin flip (Pg = ±1), and belongs to one 
of the two irreducible representations, Ag or Bu, of a 
one-dimensional lattice reflection symmetry group. 

The kinetic energy is diagonal in momentum space 
and gives rise to a cosine band, e{k) = —2tcos{k) of 
width W — At. The Coulomb repulsion is mimicked by a 
repulsive, local Hubbard interaction U, and nearest and 
next-nearest neighbor repulsions Vi, V2 . We restrict our- 
selves to the physically relevant case U > Vi > V2 > 0- 
In this work, we are not interested in issues like a com- 
plete classification of the phase diagram of the model (|^) ; 
instead, we constrain our analysis to the consideration of 
several different points in the Mott insulating phase. A 
more systematic investigation of the extended JIubbard 
model with V2 = will be published elsewheretZl. There 
it is shown that without the next-nearest neighbor inter- 
action, it is not possible to have simultaneously a small 
Mott gap and form a Mott-Hubbard exciton. 

Linear optical absorption is one of the most commonly 
used probes in experimental studies of the dynamical 
properties of a material. The optical absorption is pro- 
portional to the real part of the optical conductivity, 
which is related to the imaginary part of the current- 
current correlation function by 



II. MODEL HAMILTONIAN 



In this work psste study the one-dimensional extended 
Hubbard modcM, 



H = -t 



l-a 



Q + 1,0 



-Fi^(n,-l)(n,+i-l) 
/ 



(1) 



This Hamiltonian describes electrons with spin a =t,i 
which can hop between neighboring sites. Here c/"^, ci^^ 
are creation and annihilation operators for electrons with 
spin (7 at site Z, n/^o- = c^^ciu are the corresponding num- 
ber operators, and n; = n; ^ -f n;^|. 

Since we are interested in the Mott insulating phase, 
we exclusively consider a half-filled band where the num- 
ber of electrons N equals the number of lattice sites L. 
The lattice spacing is set to unity, ag = 1. Note that we 
have chosen the chemical potential in such a way that 
the Hamiltonian explicitly exhibits a particle-hole sym- 
metry. This Hamiltonian has two other discrete symme- 
tries which are useful for optical excitation calculations: 



(Tl(w > 0) 



Im{x(w > 0)} 



OJ 



xi^>o) = -j{o\j- — . 

L Eq - H + huj + IT] 



m 



1 



E 



l(0|ik 



tiw - {En - Eo) + IT] 



(2a) 
(2b) 

(2c) 



Here, |0) is the ground state, \n) are excited states, and 
Eq, En are their respective energies. Although 77 = 0"^ is 
infinitesimal, we may introduce a finite value to broaden 
our resonances at huj — En ^ Eq. In momentum space, 
the current operator reads 



J = 



k-(7 



sin(fc)c 



,Cfc , 



(3) 



We note that the current operator is invariant under 
the spin-flip transformation but antisymmetric under 
charge conjugation and spatial reflection. Therefore, if 
the ground state |0) belongs to the symmetry subspace 
{Ag, Pc, Ps), only excited states \n) belonging to the sym- 
metry subspace {B^, —Pc, Ps) contribute to the optical 
conductivity. According to selection rules, the matrix 
element (0|j|n) vanishes if |7i) belongs to another sym- 
metry subspace. We set h = 1 throughout, and for the 
presentation of our results we use e = i = 1 in our figures. 



2 



III. DENSITY-MATRIX RENORMALIZATION 
GROUP 

ReceBtlj^T, the density-matrix renormahzation group 
methodtSO (DMRG) has been extended. .to. the calcula- 
tion of dynamical correlation functionsOEJElJ. This nu- 
merical technique allows us to obtain iTi(tL') for all inter- 
action strengths as long as the gap is not exponentially 
small. A complete exposition of our DDMRG method 
will be published elsewherecil. 

DDMRG allows us to calculate dynamical correlation 
fimctions, such as the r.h.s. of (2b), very accurately over 
the full frequency range for fairly large systems {L < 128) 
with open boundary conditions and a finite broadening 
factor rj, i.e., the DDMRG actually gives 



S„ -Eo[u;- {E„ - EoW 



(4) 



For 77 ^ 0, ari-L{i^) reduces to ai{uj) as defined in Eq. (||). 
Ultimately, we are interested in the optical conductivity 
of an infinite system, L 00, for 77 ^ 0"*". It is shown 
in Ref. |2| that the most appropriate way of approach- 
ing this double limit is to compute a^y^L^uj) for different 
system sizes while keeping rjL = const and then to ex- 
trapolate to infinite system size. In this paper we use 
rjL = 12. 8t, which yields an energy resolution of O.lt for 
our largest system size {L = 128). 

A very useful consistency check of the method is to 
test various sum rules, relating moments of the function 
cri(w) to ground-state expectation values, which can be 
evaluated|r-with great accuracy using a standard DMRG 
methodtallS. For instance, for the Hamiltonian (|^) with 
open boundary conditions 



dio 



uja.iuj) = - {0\f\Q) 



duj 



L 

eH 
2L 



(5a) 



^(c+,Q+i,,-Hh.c.)|0), (5b) 



duJ (7i(cj) 



ni-l) 



|0) 



(5c) 



For <Tj^.i^[uj) these sum rules are not fulfilled exactly, but 
only up to errors of the order of ©([ly/t]) ov 0{[rj/t]^). 

The ground-state phase diagram of the Hamilto- 
nian (jl|) exhibits several different phases (for instance, 
Mott-Hubbard insulator, charge density wave, and bond- 
order wave for [/ > Vi > and V2 = 0, see Ref. |l^). To 
check the nature of the ground state for some fixed model 
parameters we calculate the spin gap and various ground 
state properties, such as charge density, bond order, and 
spin and density correlations, for large system sizes (up 
to L = 512 sites) with a standard DMRG method. The 
ground state of is a Mott insulator for all values of 
the model parameters used in this work. 



With DMRG one can also calculate the charge gap 
E,{L) = Eo{L + Eo{L - 1) - 2^o(i) , (6) 

where Eq {N) is the ground state energy of (Q) on a L-site 
lattice with N electrons. For L ^ 00, E^ gives the en- 
ergy threshold of the electron-hole excitation continuum. 
In a Mott insulator it corresponds to the Mott gap. In 
the one-dimensional Hubbard model (Vi = V2 = 0), it is 
known that Ec is also equal to the optical gap which we 
define as the energy threshold of the lowest band in the 
optical spectrum. In all cases with ¥±,¥2 7^ discussed 
here we have found that the optical spectrum contains a 
single band, which corresponds to unbound particle-hole 
excitations, and that E^ agrees with the onset of this 
band. Therefore, in this paper we identify the charge 
gap with the optical gap. (Some special cases for which 
the charge gap does not correspond to the optical gap 
are discussed in Ref. ^) 

We also use the symmetrized DMRG methodO to 
calculate the lowest eigenstates in the optical excita- 
tion symmetry subspace (see Sec. ||). As the standard 
DMRG method, the symmetrized DMRG yields not only 
the eigenenergies but also allows us the computation of 
various expectation values and correlation functions of 
the eigenstates (for an example, see Sec. 0). We can 
thus investigate the nature and properties of these op- 
tical excitations. In particular, it is possible to dis- 
tinguish unbound particle-hole excitations from excitons 
and from other kinds of excitations (excitonic strings, 
charge-density-wave droplets) which pan dominate the 
optical spectrum of a Mott insulatoiO. In this paper 
we consider only the regime of the Hamiltonian (0) where 
optically excited states can be described as (bound or un- 
bound) particle-hole pairs. We emphasize that the sym- 
metrized DMRG results for the optically excited states 
are always in perfect agreement with the DDMRG results 
for the optical conductivity, confirming the accuracy of 
both methods. 

All DMRG methods have a truncation error which is 
reduced by increasing the number m of retained den- 
sity matrix eigenstates (for more details, see Refs. Q^9|). 
Varying m allows one to compute physical quantities for 
different truncation errors and thus to obtain error es- 
timates on these quantities. For some quantities, espe- 
cially eigenenergies, it is possible to extrapolate the re- 
sults to the limit of vanishing truncation error and thus to 
achieve a greater accuracy. We have systematically used 
these procedures to estimate the precision of our numer- 
ical calculations and adjusted the maximal number m 
of density matrix eigenstates to reach a desired accuracy. 
The largest number of density matrix eigenstates we have 
used is TO = 1200. For all numerical results presented in 
this paper DMRG truncation errors are negligible unless 
specified explicitly. The main cause of inaccuracies are fi- 
nite size effects or extrapolation errors for L — > 00 which 
we discuss below when we present our numerical results. 
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A. Limit of large Mott gaps 

Let us now consider the situation where the Mott gap 
is much larger than the band width At. For large interac- 
tion strengths, U ^ t,Vi,V2, it is possibk-|to analyze the 
model (|l|) by means of a 1/U expansiontj. If we ignore 
corrections of the order t/U, all sites are singly occu- 
pied in the ground state. Electron transfers are limited 
to processes that conserve the number of double occu- 
pancies, and a rather simple band picture emerges for 
Vi = V2 = 0- In an optical absorption process we ex- 
cite one hole at momentum k — q/2 in the lower Hub- 
bard band, eLHB(^) = —U/2 + e(fc), and one double oc- 
cupancy at momentum k + q/2 in the upper Hubbard 
band, euHB(fc) — U/2 — e{k) (antiparallel bands). The 
total momentum of the two charge excitations is q, and 
their energy is uj. Due to spin-chsucge separation, the 
oscillator strength can be written astj 



l(Ob»|^ 



-iee{k)\ . 



(7) 



The spin sector enters the current-current correlation 
function via the momentum-dependent ground-state 
form factor gq, 



gq = 2{0\Z+^^,{q) ( ^ - S,S,+i ) ^,,,+i(g)|0) (8a) 



-{l~r) 



(8b) 



where Ts shifts all spins by one site whereas Ts> performs 
the same operation on the lattice with sites r and r + 1 
removed. 

For the large-JJ Hubbard model itself, the analysis of 
gq is rathpr involved. We can use a "no-recoil approxi- 
mation" Ej to argue that the dominant contributions to 
the conductivity come from g = and q ~ n, which cor- 
respond to vertical transitions between two antiparallel 
bands {q = 0; eLHB(fc), euHB(fc)) and between two paral- 
lel bands {q = tt; eLHB(fc), euHB(fc + 7i')). Thia-biigothesis 
has subsequently been confirmed by DDMRGollZl, which 
yields excellent agreement if the form factors are chosen 
as go = 2.65 and = 0.05 ± 0.03. Exact sum rules im- 
pose the condition 50 + = 4 ln(2) « 2.77 for an infinite 
system. The deviation of our best fits can be traced back 
to finite-size effects and numerical errors of the order of 
one percent. 

We now discuss the effects of a finite nearest-neighbor 
Coulomb repulsion Vi <C C/ with V2 = 0- It follows by 
direct inspection of the Hamiltonian (|^) that the dou- 
ble occupancy and the hole now mutually attract. To 
some extent, this is reminiscent of the situation encoun- 
tered in Wannier theory for a band insulator. How- 
ever, unlike in Wannier theory, the double occupancy 
and the hole are not fermionic quasi-particles but spin- 
less hard-core bosons. Even more importantly, there is a 
critical value Vc — 2t below which no exciton appears 



below the threshold of the particle-hole continuum at 
ijjph = E(. — U — At. For Vi > Vc there is an exciton 
at the energy 

= C/ - 14 - At^/Vi . (9a) 

In addition, there is a second Mott-Hubbard exciton at 
the energy 



W2 = [/ - Vl , 



(9b) 



which, for Vi < 4i, lies in the particle-hole continuum 
but carries very little spectral weight. 
The optical conductivity is given by 

uj(Ji{uj) — ■Kg-,^t'^e^5{uj — 002) 

+gQt^e^\Q{Vi - 2t)n (l - {2t/Vif) S{u; - ui) 



2+2 /T 

+Q{At ~ \io ^ U\) ^ 



iiu^'U)/ity 



Vi{uj - ui) 



(10) 



Here, Q{x) is the Heaviside step function. Apart from 
the two (5-peaks corresponding to the excitons, there is 
a particle-hole continuum for jcj — [/| < At. Near the 
lower (upper) boundary, the optical conductivity shows 
a characteristic square-root increase (decrease). The only 
exception is the case of Vi = Vc where the optical con- 
ductivity diverges at the threshold, 



(Tl(ljj) C>C 



for UJ — t 



0^ 



(Vi = 2t) 



(11) 



The predictions of this strong-coupling analysis are 
confirmed by our numerical results. Note that ai(uj) ^ 
1/U so that it is more convenient in DDMRG to calculate 
directly the imaginary part of the current-current corre- 
lation function, i.e., Im{x(w)} = wcri(w). In analogy to 
Eq. (Q), a broadening 77 is introduced in the DDMRG 
procedure for the current-current correlation function, 



Im{x,,;L(t^)} = 7X1 



[lu - [En - EoW 



{12) 



For 7y — > 0, this expression reduces to Im{x(a;)}, and 
we analyze the finite-size effects as discussed above for 

(Tl(w). 

Figure |l| shows Im{x^_L('^)} for L = 128 and rj = Q.lt 
obtained in the large-?/ limit of the extended Hubbard 
model with Vi = 5t and V2 — 0. We compare the 
DDMRG data to the analytical formula (10), convolved 
with a Lorentzian of width ry. The DDMRG data and 
the large-C/ result are in very good agreement when we 
choose go — 2.65, = 0.08, as discussed above. These 
values are found to be essentially independent of Vi . We 
previously obtained a similarly good agreement far the 
case of the laige-U Hubbard model {Vi = V2 ~ 0)l3. 
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FIG. 1. Current-current correlation function lm{x{i^)} for 
U/t ^ oo, Vi = 5t, Vi = 0, and 77 = O.lf. The solid line is 
the DDMRG result Im {x,,;L(t^)} for L = 128. The dotted 
line is Eq. ( p^ convolved with a Lorentzian of width r\. Note 
the logarithmic scale of the ordinate. Inset; same results on 
a linear scale. 

For V\ = 5t, V2 = 0, most of the spectral weight is 
carried by the exciton atw — [/ — uJi — U — —5.8t, as one 
also observes for typical band insulators. Therefore we 
use a logarithmic scale to make visible the contributions 
of the other exciton and the particle-hole continuum. The 
use of a logarithmic scale also reveals deviations around 
the second exciton {uj — U ^2 — U = —5t) and close 
to the upper band edge (uj — U fa +4t), which are as- 
sociated with difficulties in the numerical determination 
of small contributions to the optical spectrum, and re- 
maining finite-size and boundary effects. In fact, these 
deviations are less than one percent of the total spectral 
weight, and are completely irrelevant for practical, i.e., 
experimental purposes. On a linear scale they are not 
visible as seen in the inset of Fig. |l|. 

The case of < V2 < Vi can be treated analogously 
and does not yield any new physical aspects. It is known 
that increasing the interaction range simply reduces the 
critical coupling below which no exciton appears in the 
gap and pjcreases the number of visible excitons in the 
spectruma. 



B. Regime of intermediate Mott gaps 

In the simple Hubbard model {Vi = V2 ~ 0) we 
previously found that the optical conductivity evolved 
smoothly from the regime of small Mott gaas {U ^ t) 
to the limit of large Mott gaps {U ^ t)cj. Optical 
excitations can simply be interpreted as a particle-hole 
pair, i.e., a pair of spinless quasi-particles with opposite 
charges, corresponding to the hole and double occupancy 
for large Mott gaps and to the holon-antiholon pair in the 
field-theory limit of small Mott gaps (see Sec. IV). 

The optical spectrum of the extended Hubbard 
model (|l|) shows more diversity in the presence of a finite 
intersite Coulomb repulsion (Vi > V2 > 0). Even within 



the Mott insulator phase one can observe "exotic" optical 
excitations S|WWj;h as charge-density-wave droplets or exci- 
tonic stringstD. In this study we restrict ourselves to the 
Mott insulator regime of (|]) where the dominant optical 
excitations can be described (bound or unbound) 

particle-hole pair and the Coulomb interaction is strong 
enough to generate at least one bound pair (exciton). 

Varying the model parameters U/t, Vi/t, and V2/t we 
have investigated the optical excitations of systems with 
a Mott gap ranging from 10 to 0.1 times the band width 
4t. It is important to note that the Mott gap Ec incrpaps 
with U and V2 but decreases with increasing ViOJ'EZlEJ. 
As in the large Mott gap limit, we have found that the 
intersite Coulomb interaction must exceed a critical value 
before a discrete absorption peak appears at an energy 
Wex below the optical gap Ec- For V2 = the critical 
value is Vi « 2t for all U/t in agreement with ourjapalyt- 
ical strong-coupling analysis and a previous workcJ. The 
critical value oiVi/t becomes smaller as the next-nearest 
neighbor repulsion V2 increases. 

We have analyzed the nature of the excited states as- 
sociated with the discrete absorption peak using various 
measurements and correlation functions. For instance, in 



Sec. VB we present a method to determine the size of a 
particle-hole pair. This analysis confirms that this ex- 
cited state is clearly a bound particle- hole pair (exciton). 
The exciton binding energies SE = Ec — observed 
in our calculations range from SE — Q.03t to SE = 12t 
and the exciton sizes measured with the procedure of 
Sec. VB vary from 20ao down to slightly more than one 



lattice spacing. 

It is interesting to note that we have never found 
more than one exciton in the regime of intermediate 
Mott gaps. Our numerical results (DDMRG and sym- 
metrized DMRG) for finite open chains sometimes yield 
more than one optical excitations with energy coq ^ lucx + 



c{L)q < Ec and quasi-momc 



0,1,2, 



1 ^■n 



jta q w 7r(2£-h 1)/(L + 1), 
3 The first of these states 



{£ = 0) has always much more spectral weight than the 
other ones and corresponds to an exciton with momen- 
tum 5 = in an infinite chain. Scattering by the chain 
ends is responsible for the small but finite optical weight 
of the other states {£ > 1), corresponding to excitons with 
momentum (7 7^ in an infinite system. Thus, with pe- 
riodic boundary conditions or in an infinite system, only 
one exciton contributes to the optical conductivity ai (w) 
of the Hamiltonian (|l]) in the regime of intermediate Mott 
gaps. 

In contrast to this, both the strong-coupling analy- 
sis and field theory allow for more than one exciton in 
the optical spectrum of a Mott insulator in the thermo- 
dynamic limit if the Coulomb repulsion becomes strong 
enough. For the model (Q) an increase of Vi and V2 does 
not lead to the formation of a second Mott-Hubbard 
exciton. Instead, the nature of the lowest optical ex- 
citations changes to charge-density-wave droplets or ex- 
citonic strings, or the ground state develops long-range 
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order. Both the strong-couphng analysis and the field- 
theory approach assume a Mott insulator ground state 
and particle-hole pairs as optical excitations, and thus 
do not reproduce this instability toward charge density 
ordering. It is conceivable, though, that the inclusion 
of Coulomb terms beyond next-nearest neighbors in the 
lattice model favors the appearance of more excitons 
in the optical conductivity of a Mott insulator. 

Besides the single exciton peak we have always found 
that ai{u}) shows an absorption band starting at the 
charge gap E^- Within the resolution of our method the 
optical spectrum does not display any other feature. The 
investigation of the excited states in the continuum above 
Ec is much more demanding than the analysis of the iso- 
lated exciton peak. Whenever this has been feasible, we 
have found that the excited states contributing to the 
absorption continuum can be described as an unbound 
particle- hole pair (see Sec. |VB| ). 

As a typical example, the optical conductivity of (|l|) 
for U = 8t, Vi = 4t, and V2 = 2t is shown in Fig. |. A 
peak at the exciton energy Wex = 3.34t is the dominant 
feature of the spectrum while a very weak band is visi- 
ble for uj > Ec = 4.05t. No gap is visible between Wox 
and Ec because of the broadening of the strong exciton 
peak. In the inset of Fig. || one can see the weak particle- 
hole continuum part of the spectrum separated from the 
strong exciton contribution. The onset of the absorption 
band is clearly visible at a; « = 4.05t. The small ir- 
regular fluctuations seen in the inset are numerical errors 
(truncation errors) made visible by the small scale used. 

In summary, our numerical results show that there is 
no qualitative change in the optical conductivity of a 
Mott insulator with excitons when one goes from the 
limit of large Mott gaps down to the regime of inter- 
mediate Mott gaps with Ec ^ QAt. As in the Hubbard 
modelly, the simple picture of the strong-coupling anal- 
ysis remains qualitatively valid. 




CO 

FIG. 2. Optical conductivity Gr,-L{sjj) for U = 8t, Vi = 4f, 
and V2 = "it, calculated with DDMRG on a 128-site lattice 
(77 = Q.lt). The inset shows an expanded view of (t,,;l(i^) 
(dotted line) for 3.5 < oj/t < 6. The exciton (dashed) and 
continuum (solid) contributions to cr^;L(cij) are also shown. 



IV. FIELD THEORY 

In order to address the regime of small Mott gaps, 
we study the extended Hubbard model (^ in the field- 
theoretical limit. This limit can be constructed di- 
rectly from the lattice model in the weak-coupling regime 
U,Vi,V2 ^ t. The low-energy physics of the non- 
interacting model is simply described in terms of a mass- 
less Dirac fermion with velocity vp = 2taQ. The inter- 
actions introduce a four-fermion coupling. The resulting 
effective theory is kno.wn as the U{1) Thirring model and 
can be represented asEZrnJ 

[ dx[nc + Hs] , 

Uc 



(13) 



where J", J" (/", /") are left and right moving SU(2) 
spin currents (SU(2) pseudospin currents) and 

gi_ = 2{U - 2[Vi - V2])ao , g\\ = 2(C/ + 6^1 + 2^2)^0 , 



, 2iU/A + Vi+V2)ao 

VF H , Ws 



VF 



The Hamiltonian (|T^) explicitly exhibits spin-charge sep- 
aration: 7ic,s describe the charge and spin degrees of free- 
dom, respectively, which are independent of one another. 
As long as g± > 0, the spin sector remains gapless and 
can be bosonized in terms of a Gaussian model. The 
charge sector can be bosonized as well, as is for exam- 
ple shown in Ref. |2^. The result is a Sine- Gordon model 
(SGM) 



lOTT 



(vcdx'pc)'^] + 2^cCos/ 



(15) 



where fic and (3 are functions of U, Vi, ¥2- In order to 
utilize results obtained from the integrability of the Sine- 
Gordon model we use a flat renormalization scheme, in 
which /3 is a constant 



Attvf 



Attvf 



9t\ 



9l 



(16) 



The pure Hubbard model corresponds to the limit /3 — s- 1 
and the effect of Vi and V2 is to decrease the value of (3. 

In the field-theory limit the electrical current operator 
is found to be 



J{t,x) = ^/A'^tC^c 



(17) 



where A' > is a non-universal constant. As seen from 
Eq. (|l^) , the current operator does not couple to the spin 
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sector. This shows that spinons do not contribute to the 
optical conductivity. Therefore we ignore the spin part 
Tis of the Hamiltonian from now on. 

The calculation of the optical conductivity in the field- 
theory limit has thus been reduced to the evaluation of 
the retarded current-current correlation function in the 
charge sector, 

X^^M = — / dtexp{iiut) dx{[J{x,t),J{0,0)]). 



LaQ Jq 



(18) 



We turn to the calculation of this correlation function in 



Sec. [VA 



For the-, Sine-Gordon model, many exact results are 
availablecj. The spectrum of the SGM depends on the 
value of the coupling constant (3. In the so-called repul- 
sive regime, 1/^/2 < /? < 1, there are only two elemen- 
tary excitations, called soliton and antisoliton. From the 
point of view of the underlying lattice model (|l|) these 
correspond to holon and antiholon (spinless excitations 
of opposite charges). These have a massive relativistic 
dispersion. 



E{P) = y^lVP + V^P^ , 



(19) 



where M is the single-particle gap which is related to the 
optical gap by A = 2M. At weak coupling the gap scales 
as 



M 



8t 



igx+2}/4gx 



(20) 



where we have fixed the constant factor by comparing to 
the exact result for the Hubbard model, and where 



U + 6V2 + 2V2 



{U + Wl+2V2)l2TTt . 



1/2 



(21) 



We note that the gap vanishes on the critical surface 
U — 2Vi + 2V2 = separating the Mott-insulating phase 
with gapless spin excitations from another phase with a 
spin gap. 

In the regime < /? < l/\/2, soliton and antisoliton 
attract and can form bound states. In the SGM these are 
usually known as "breathers" and correspond to excitons 
in our lattice model. There are 



1-/3' 



21 



(22) 



different types of excitons, where [x] in ( p2| ) denotes the 
integer part of x. The exciton gaps are given by 



M„ = 2M sin(n</2) , n = 1, , 



(23) 



1-/32 



(24) 



Therefore, the single-particle gap M and the coupling f3 
fully characterize the spectrum of the SGM. 

One knows that the field-theory approximation to the 
lattice problem is valid in the limit U,Vi,V2 <C t where 
the single-particle gap is much smaller than the band- 
width—For the Hubbard model {Vi — V2 — 0) we have 
foundEfl that field theory gives surprisingly good results 
for the optical conductivity even for intermediate single- 
particle gaps of magnitude M « 0.3t. The same holds 
true for the Mott insulating phasii.p£ the extended Hub- 
bard model with Vi > 0, V2 ^ OEJEJ. As we shall show 



in Sect. tVB, field theory remains applicable even in the 



presence of excitons. 

In the framework of the field-theory approximation to 
the lattice problem we can determine the value of P only 
in the hmit U, Vi, P2 ^ In fact, for the-Hubbard model 
/3 = 1 is fixed by the S0(4) symmetryEd. As seen from 
Eq. (p^, the effect of a small Vi, V2 is to decrease the 
value of p. One may therefore hope that by carefully 
tuning U, Vi and V2 one may stay in a regime with a 
"small" single-particle gap, i.e., close to the critical sur- 
face U — 2Vi + 2V2 = 0, but with a sufficiently small P for 
excitons to exist. We have found numerically that it is 
indeed possible to reach regions of the parameter space 
where field theory is valid and P'^ is as small as 0.36 close 
to a critical surface U — 2Vi + 2V2 ~ which separates the 
Mott insulating phase from other phases with long-range 
order. The determination of the field-theoretical parame- 
ter /J^ as a function of the lattice model param eters using 
DMRG results will be discussed in Sec. IV B. 



where 



A. Optical conductivity in the Sine-Gordon model 



Our task is now to calculate the Fourier transform of 
the retarded, dynamical current-current correlation func- 
tion (|l|) in the Sine-Gordon model. This is done by go- 
ing to the spectral representation and then utilizing the 
integrability of the SGM to determine exactly the ma- 
trix elements of the current operator between the ground 
state and various excited states. This piettod is known 
as the form-factor bootstrap approachcSH and has re- 
cently been applied to calculate th e opt ical conductivity 
in the repulsive regime of the SGMllal23. Here we review 
briefly the relevant steps and refer to Ref. ^ for a more 
detailed exposition. In order to utilize the spectral rep- 
resentation we need to specify a basis of eigenstates of 
the Hamiltonian (15). Such a basis is given by scattering 
states of excitons, holons and antiholons. In order to dis- 
tinguish these we introduce labels ei, 62, . . . cat^^, h, h. As 
usual for particles with relativistic dispersion, it is useful 
to introduce a rapidity variable 9 to parameterize energy 
and momentum 
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Eh{e) = M cosh 9 , Ph{0) = {M/vc)smh6 , (25a) 
Ef,{9) = M cosh 61 , Pf,{e) = {M/vc)smhe , (25b) 
{9) = Af„ cosh e , Pe„ (8) ^ {Mn/vc) sinh 6 , (25c) 



where the exciton gaps M„ are given by (|2^). Next we 
turn to the construction of a basis of scattering states of 
holons, antiholons and excitons. A convenient formahsm 
to this end is is obtained by using the Zamolodchikov- 
Faddeev (ZF) algebra. The ZF algebra can be consid- 
ered to be the extension of the algebra of creation and 
annihilation operators for free fermions or bosons to the 
case or interacting particles with factorizable scattering. 
The ZF algebra is based on the knowledge of-the exact 
spectrum and scattering matrix of the modeo. For the 
SGM the ZF operators (and their hermitian conjugates) 
satisfy the following algebra 



(26a) 



zli02)s:f;i02~0i)z^Hei) 

(27r)C,S(0i-02) . 



:,e2,<!i. 



(26c) 



Here Sl}'l?{0) are the known (factorizable) two-particle 

scattering matrices0 and Sj — h,h, ei, . . . , e[i/^]. 

Using the ZF operators a Fock space of states can be 
constructed as follows. The vacuum is defined by 



(27) 



Multiparticle states are then obtained by acting with 
strings of creation operators Z^{0) on the vacuum 



\0n ■■ - 91) ^ zlie,,) ... zlie^m. 



(28) 



In term of this basis the resolution of the identity is given 
by 



1 = 

CX 

+EE 

n— 1 ei 



(29) 



d0i... d0„ 
(27r)"n! 



h)e„...e,''-'"{9l...0n\ 



Inserting (g^ between the current operators in (^, we 
obtain the following spectral representation of the corre- 
lation function 



d 0i... d0n 

{2n)"nl 

n 

exp(i^ - E,t^\{0\J{0, O)\0n . . . 0i)e^...e, W (30) 



where Pj and Ej are given by 



Pj = — smh0j, Ej = cosh6'j , 



(31) 



and 



/ 



(01^(0,0) |0„., 



n/e„...ei 



(32) 



are the form factors. Our conventions in (31) are such 
that Mi^ — Mj^ — M and M,,^ — Af„. After carrying out 
the double Fourier transform we arrive at 



-27r 
Lao 



EE 



n— 1 ei 

5{T.jM,^smh0j/v, 



d0i . . . d0n I J. J 

1/ 



(27r)"n! 



ijj - X] o M^. cosh 0.J +iT] Lu + J2j cosh 0j + irj 



(33) 



This then yields the following representation for the real 
part of the optical conductivity {uj > 0) 



Zl [0,)Zl {02} = Zl {02)Zl [0,)Sli:ll [0, - 02) , (26b) ^"^^ 



2^ 

LaQLO 



EE 



71—1 ei 



X |/-^(ei 



i)ei. 



d0i... d0„ 

(27r)"n! 

|2 



(34) 



M 



S{Y^ smh0k)S{uj - Me, cosh6'fc) . 



The missing ingredient in ( p4[ ) are the form factors. In 
Refs. |3^j3^ integral representations for the form factors 
of the current operator in the Sine-Gordon model were 
derived. Using these results we can determine the first 
few terms of the expansion (^). In particular, the form 
factors involving excitons are determined via tiie boot- 
strap axioms for soliton-antisoliton form factorg23. 

From (^^ it is easy to see for any given frequency 
Lu only a finite number of intermediate states will con- 
tribute: the delta function forces the sum of single- 
particle gaps J2j -^(j to be less than cj. Expansions of 
the form ( |3^ ) are usually found to exhibit a rapid conver- 
gence, which-i^n be understood in terms of phase space 
argumenta^'Lj. Therefore we expect that summing the 
first few terms in ( p^ ) will give us good results over a 
large frequency range. 

Using the transformation property of the current op- 
erator under charge conjugation one finds that many of 
the form factors in ( |3^ ) actually vanish. In particular, 
only the "odd" excitons ei, 63, . . . (assuming they exists, 
i.e., /3 is sufficiently small) couple to the current oper- 
ator. The first few non-vanishing terms of the spectral 
representation (p4) are 



^FT 



(c) 



A 



[(Wex + l)/2] 

E ^e2„-i('^) +cr/,l(^) +CTeie2(^ 



(35) 



Here A = A'vc/Lao and (Te„(u;), cr^^(w) and tTeie2('^) are 
the contributions of the odd excitons, the holon-antiholon 
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continuum and the 6162 exciton-exciton continuum re- 
spectively. The latter of course exists only if N^x > 2. 
We find 



ML-1 



f2n-l5{uj - M2„_l) 



(36a) 



(36b) 



n=l 



X exp —2 



dt sinh(t(l - $)/2) sinh^(mt^/2) 



t sinh(<^/2)cosh(t/2) 



sinht 



The holorbiantiholon contribution has previously been 
determinedEEl and is given by 



X exp 



4Vw2 - 4M2e(w ~ 2M) 

u;2[cos(7r/C) +cosh(e'/0] 

dt sinh[i(l - ^)/2] [1 - cos{t9/Tr) cosht] 



(37) 



sinh(i^/2) cosh(t/2) sinht 



where 9 — 2arccosh(ti;/2M). Finally, we quote the result 
for the 6162 exciton-exciton continuum 



2co\h2{0l2)\' 



Ml 



Mir 



{2MIM2) 



(38) 



where 

|/l2(^?12)|' =A' 



2 ^gitauTT^I sinh^ 012 + sin^(7r^/2) 
"2 sinh^ 6(12 +sin2(37rC/2) 



ex ( 8 f°°'^* ^^"^ ^ smh{t£_) sinh[t(l ^1 cosh(2if ) \ 
^""^V X ~ sinh2(2t) j 

' sinh(2t^) sinh[i(l + ^)] cos(2i6li2/C) 
t coshf sinh(2/;) 



X exp |— 4 

X (4cos [27r(cosh6'i2 + cos«/2))/^])" 
and 



(39) 




For 1 > > 1/3 the contributions of the first odd 
exciton ( ^ ) and of the holon-antiholon continuum (37) 
dominate the field-theoretical optical conductivity (35) 
(at least in the low-frequency regime) and other contri- 
butions vanish or can be neglected. The optical conduc- 
tivity can be written as 



(41) 



where hj3{x) is a universal function depending only on 0^ . 
Therefore, only the parameter 0^ determines the shape 
of the optical spectrum. The parameters A and A/ A just 
set the energy and conductivity scales, respectively. 



For 1 > > 1/2 the optical spectrum contains a sin- 
gle band, while for 1/2 > 0^ > 1/3 it contains one band 
and one exciton peak in the optical gap at the energy 
LUcx — Ml . The optical weight is progressively transfered 
from the band to the single exciton as decreases down 
to 1/3. The absorption band increases smoothly at the 
threshold A, as 



^FT 



(lu) oc \/w — a for u! ■ 



(42) 



for all values of (3'^ but /J^ — 1/2, when the exciton 
peak merges with the band. In this case af^lui) shows a 
square-root divergence at the absorption threshold 



ai (uj) (X 



1 



OJ - 



for Lo^ A+ {13^ = 1/2) . (43) 



The behavior of the field-theoretical conductivity is thus 
qualitatively similar to that of the lajttice model in the 
large Mott-gap limit (|l^). For V2 = dlZl it has also been 
found for intermediate Mott gaps that cti(ijj) diverges as 
a square-root at the absorption band threshold for the 
critical coupling V = Vc below which no exciton appears 
in the optical gap. For Vi 7^ Vc, cri{uj) increases smoothly 
at the absorption band threshold. This generic behavior 
of cri(uj) illustrates once more the absence of significant 
qualitative changes in the particle-hole continuum and 
single-exciton spectrum of (|l|) as one goes from the large 
to the small Mott gap hmit. For < 1/3, the field- 
theoretical optical conductivity shows more features, but 
it seems that this regime cannot be reached in the lattice 
model (!]) and thus we shall not discuss it further. 



B. Application to the lattice model 

The field theory optical conductivity af'^{uj), (Eq. ^5| ), 
depends on three parameters: /J^, the gap A = 2M, and 
the normalization A. To compare the field theory predic- 
tions with our numerical results for the lattice model (|l|) 
one needs to determine the field theory parameters cor- 
responding to specific values of the model parameters U, 
Vi , and V2 . Fortunately, if there is exactly one exciton in 
the spectrum (1/2 > > 1/3), this can be done using 
standard DMRG techniques. 

The first step is the calculation of the gap parame- 
ter A. The gap A obviously corresponds to the charge 
gap (^ extrapolated to the infinite system limit, A = 
\m\L^oo Ec{L). The second step is determining /J^. The 
exciton energy ujc-x{L) of the lattice system can be cal- 
culated using the symmetrized DMRG method. The ex- 
trapolation to infinite system size gives us the exciton gap 
of field theory, Eq. (||), Mi = limL^oo Wcx(i). The pa- 
rameter P'^ is then fixed by the ratio of Mi and A = 2M 
through Eqs. and (||). 

The final step is the calculation of the normalization A. 
Evaluating A accurately turns out to be the most difficult 
task and we have tried two different approaches. In the 
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first one we calculate the excitonic spectral weight in the 
lattice model using the symmetrized DMRG technique 
and extrapolate to an infinite system size. This yields 
the normalization A by comparison with the field-theory 
prediction, Eqs. ( |35| ) and (^6|). This method is simple 
and exact on the field theory side but difficult to carry 
out numerically because of significant truncation errors 
and complicated finite-size effects. 

In the second approach we use the sum rules (|). We 
numerically integrate the contributions of the exciton and 
the holon-antiholon continuum to the optical conductiv- 
ity for a fixed value of The result has a triv- 
ial dependence on A and A because of (41). This gives 
us the l.h.s. of (|5bD or (g) assuming that the neglected 



contributions toai"^{uj) are insignificant and that most 
of the optical spectrum weight is concentrated at low en- 
ergy where field theory is expected to describe the lattice- 
model properties accurately. The r.h.s. of ( |5b| ) and (^) 
can be calculated for the lattice model with DMRG and 
extrapolated to the infinite system limit. Comparison of 
both sides of the equations give the value of the normal- 
ization A. 

We prefer the second approach because it is simpler 
and more accurate than the first one as far as the numer- 
ical calculations are concerned. It can also be used when 
there is no exciton (1 > > 1/2). Unfortunately, the 
second approach works only if the exciton (|3^) and holon- 
antiholon continuum ( ^7|) reproduce the lattice model 
ai{uj) very accurately, i.e., if 



SSn 



dujuj-'^ \Aae^{uj) + Aaf^j^Lo) - (Ji{uj)\ (44) 



is very small compared to the l.h.s of ( pb| ) for n = or 
the l.h.s. of ( |5c| ) for n = 1. 

The validity of the field-theory approximation to the 
model (|l|) is not necessarily restricted to the limit 
U,Vi,V2 <C t but rather to the regions where the single- 
particle gap M — A/2 = Ec/2 is smalLcpmpared to the 
band width U. In the Hubbard modelliJ {Vi = V2 = 0) 
we have found that field theory describes the optical con- 
ductivity accurately even for U = 3t, corresponding to 
A fa 0.6t. In the case Vi > 0, V2 = 0, field theory is valid 
not only in the weak-coupling limit {U, V\ ^ t) but also 
in the vicinity of a critical line U — 2Vi ~ separating 
the Mgtt-insulating phase from phases with long-range 
ordertZl. The Mott gap vanishes or becomes extremely 
small on this critical line at least up to [/ = At. The 
holon-antiholon continuum contribution (^ ) to the op- 
tical conductivity agrees with DDMRG results for gaps 
A « 0.6t and ranging from 1 to 1/2. 

In the general case (?7 > Vi > V2 > 0), there is also a 
critical surface U — 2Vi + 2V2 ~ separating the Mott in- 
sulating phase from other phases (in particular, a charge- 
density- wave phase), where the Mott gap vanishes or be- 
comes very small even for relatively strong Coulomb in- 
teraction (at least up to C/ = 6t). A careful tuning of 
the parameters U, Vi, and V2 allows one to reach re- 
gions of the parameter space, where field theory appears 



to remain valid and decreases down to 0.36 (for in- 
stance, U ~ 6t, Vi — 4.5t, and V2 = 2t). We can thus 
compare our numerical results with the field-theoretical 
predictions for the optical conductivity in the parameter 
regime with a single exciton. Again, we have found a 
good agreement between DDMRG and field theory for 
gaps as large as A = 0.6t. 

In Fig. ^ we compare the optical conductivity from 
field theory and DDMRG for U = 5.2t, Vi = 3.7t, and 
V2 = 1.3t. Only the low-frequency {ui < 2t) part of the 
spectrum is shown as there is almost no spectral weight 
at higher frequency. The numerical results have been cal- 
culated on a 128-site chain with a resolution of 77 = OAt. 
The field-theory result, Eqs. (|35|)-(p7|), has been con- 
volved with a Lorentzian of width rj = OAt to allow for a 
direct comparison with DDMRG results. The field theory 
parameters A = 0.625^, = 0.40, and A = lA2e'^t have 
been evaluated for an infinite chain as discussed above. 
The gap between the exciton peak at Mi = 0.544i and 
the threshold of a weak holon-antiholon continuum at A 
is not visible in Fig. |^ because it is smaller than the finite 
resolution introduced by the broadening. 

In Fig. H one sees that the agreement between numer- 
ical and field-theory results is good. The visible discrep- 
ancies are well-understood finite-size effects: the exci- 
tation energies of (^), in particular the exciton energy, 
decreases as the system size increases, the total spectral 
weight is slightly smaller in the finite chain as shown by 
corrections to the kinetic energy (r.h.s. of (|5b|)) of the or- 
der 1/L, and the exciton peak is broadened and flattened 
because of the scattering by chain ends. If we choose the 
field-theory parameters A and A (the energy and conduc- 
tivity scales) to fit the finite-system DDMRG conductiv- 
ity, differences almost completely vanish for A — 0.669i 
and A = 1.02e^i, as seen in the inset of Fig. ||. 




FIG. 3. Optical conductivity cr^;L(aj) (solid) calculated 
witli DDMRG for U = 5.2t, Vi = 3.7t, and V2 = 1.3t on 
a 128-site lattice {rj — O.lt). Field-theoretical result aj"^ (lo) 
(dashed) for A = 0.625i, = 0.40 and A = l.Ue'^t, con- 
volved with a Lorentzian of width rj — O.lt. Inset: same 
results with A = 0.669t and A = 1.02e^t. 
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As mentioned in Sec. IIIB we have never observed 



more than one exciton below the optical absorption con- 
tinuum in the regime of intermediate Mott gaps, down to 
Ec = OAt. Therefore, we cannot evaluate the field-theory 
predictions for the additional excitons and the exciton- 
exciton continuum in the regime < I. It remains 
conceivable, however, that smaller values of can be 
reached in the model (0) in the limit = A <^ t. 



V. EXCITON PROPERTIES 



In Wannier theoryu exciton properties, such as size, 
binding energy, effective mass, or optical weight, are re- 
lated by simple equations and exhibit a monotonic behav- 
ior as a function of the Coulomb repulsion strength. This 
simplicity is due to a drastic assumption of this theory: 
optical excitations are represented by an electron (in the 
conduction band) and a hole (in the valence band), which 
are completely independent from the system's other de- 
grees of freedom. The interaction with these degrees of 
freedom is taken into account only through renormalized 
parameters such as effective masses for the electron and 
hole, and an effective background dielectric constant for 
the Coulomb interaction between electron and hole. 

In a Mott insulator the exciton properties show a more 
complex behavior. Although a Mott-Hubbard exciton 
can also be described as a bound pair of excitations with 
opposite charges, the Coulomb interaction at the same 
time determines the size of the Mott gap, the exciton 
properties, and the couphng of the exciton to the other 
electrons in the system. Therefore, an increase of the 
Coulomb interaction strength does not simply bind the 
exciton more tightly, but also renormalizes the gap and 
couples the particle-hole excitation more strongly to the 
other electrons. This leads to a non-monotonic behav- 
ior of Mott-Hubbard exciton properties as a function of 
the Coulomb repulsion strength, and even to an instabil- 
ity toward the formation of excitonic strings or charge- 
density- wave dropletsEZl. 

In the following two subsections we shall discuss the 
binding energy and size of a Mott-Hubbard exciton us- 
ing analytical results in the limits of large and small Mott 
gaps, and numerical results in the intermediate regime. 



A. Binding energy 

The exciton binding energy is usually defined as the 
energy difference between the exciton Wex and the band 
edge of the particlc-holc continuum Ec 



SE = Er-< 



(45) 



LU 

to 




m 



- 1 



In Wannier theory this quantity is the energy required to 
break an exciton into independent hole and electron. 



FIG. 4. Charge gap Ec (circles, right axis) and exciton 
binding energy 5E (squares, left axis) as a function of Vi 
for U = Gt and V2 = 2t. 

In a Mott insulator this is not always warranted as 
6E does not necessarily correspond to the minimal en- 
ergy required to break an exciton. For instance, as seen 
in (9b), there is an exciton above the band edge in the 
strong-coupling limit U ^ At > Vi (V2 = 0). Except for 
this strong-coupling limit, however, we have only found 
excitons with an energy cJqx lower than the absorption 
band edge Ec in our study of the model (^. Therefore, 
we assume that 6E is indeed the binding energy of the 
single Mott-Hubbard exciton observed in the spectrum 
of the Mott-insulating phase of (Q). 

In the limit of a large Mott gap {U ^ t,Vi,V2) 
5E ^Vi+ At^/Vi - it for Vi > 2t and V2 = 0, see (|9|). 
The binding energy increases monotonically with Vi but, 
clearly, it is not a good measure for the strength of the 
Coulomb interaction as it is independent of U and van- 
ishes even when Vi is still quite strong (Vi = 2t). For in- 
termediate Mott gaps (0.1 < Ec/At < 10) we have found 
that the binding energies (ranging from 0.03t to 12t) do 
not depend in a simple way on the Coulomb interaction 
strength. In particular, SE has a non-monotonic behav- 
ior as illustrated in Fig. ^. The binding energy 6E first 
increases with Vi as long as the gap Ec remains essen- 
tially unchanged, then SE and Ec decrease rapidly as 
Vi approaches the critical surface U — 2Vi + 2Vi ~ 0. 
It is likely that the same behavior also holds for small 
Mott gaps {Ec — A t), where field theory gives 
5E = A(l - sin(7r^/2)), ^ < 1, for the first exciton, 
see (H^). If we take the weak-coupling result (|l^) as an 
indication of the qualitative dependence of /? on Vi , then 
an increase of Vi leads to smaller ^, (|2^), and thus to a 
larger binding energy. On the other hand, an increase of 
Vi can also leads to a sharp decrease of the gap A close 
to the critical surface U — 2Vi + 2V2 « 0, if we again take 
the weak-coupling results (|2^) and (|2^) to be indicative 
of the qualitative dependence of A on • Depending on 
the other parameter values, one of these two effects dom- 
inates and leads to an increase or decrease of the binding 
energy when Vi becomes larger. 

In summary, our analysis shows that the binding en- 
ergy of a Mott-Hubbard exciton does not provide a good 
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estimate of the Coulomb interaction strength in a Mott 
insulator. Of course, a large gap Ec or binding energy SE 
requires a strong Coulomb interaction. In general, how- 
ever, a small Ec or SE are not an indication for a weak 
electron-electron interaction. In contrast to the views of 
Rcf. even a "small" exciton binding energy, of the or- 
der of 0.1 eV in some 7r-conjugated systems, does not im- 
ply that electron-electron interactions are small in these 
materials. 



B. Size 

In the limit of large Mott gaps {U ^ t,Vi,V2) an 
optical excitation is simply made of a hole and a dou- 
ble o ccupa ncy in a background of singly occupied sites 
(Sec. [II A ). The probability of finding the hole and dou- 
ble occupancy at a distance x in an optical excitation is 
given by the correlation function 



Chd(a;) = (ni^i;ni^i{l-hi+^,^){l~ni+^^i)) , (46) 

where (...) means the expectation value in the corre- 
sponding excited iV-electron eigenstate. The average dis- 
tance between hole and double occupancy is 



Chd 



j:.cux)\x\ 



(47) 



If the hole and double occupancy are not bound, this 
quantity diverges as the system size L goes to infinity. 
For an exciton, Chd remains finite as L —^ oo and can be 
interpreted as the exciton size. For Vi > 2t and V2 = 
an analytical calculation gives the exact result 



Chd{x ^ 0) cc exp(— k|2:|) 



(48) 



with K = 21n(Vi/2i) for the lowest exciton. The exciton 
size is then 



Chd — 



1 



1 



1 - {2t/viy 



(49a) 
(49b) 



Chd diverges as Vi approaches the critical value (I4 = 2t), 
below which the pair is not bound, and tends to unity 
for strong coupling, pij^2i 00. Using a density-density 
correlation functionBll3 yields equivalent results. 

Unfortunately, the correlation function for hole and 
double occupancy and the density-density correlation 
functions yield unclear results in the regime of intermedi- 
ate Mott gaps because the ground state already contains 
a finite density of holes and double occupancies when 
Ec/t is finite. The function (^), for instance, tends to 
a finite value cP as a; — + 00, which is determined by the 
density of doubly occupied sites d = {fii^-^fii^i). These 
quantum charge fiuctuations hide the weak exciton con- 
tribution to (E|). 
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FIG. 5. Correlation functions for a hole and double occu- 
pancy, Chd (a;) (dashed), and for an electron- hole excitation, 
Cch{x) (solid), calculated for the lowest exciton in a 96-site 
system with U = 40t, Vi = 2.5t, and V2 = 0. The center of 
the pair is in the middle of the lattice. 

In Fig. ||we show the correlation function C\^d{x) cal- 
culated with DMRG for the lowest optically excited state 
in the model (g) with U = 40t, Vi = 2M, and V2 = 0. 
This system is in the regime of large Mott gaps with 
Ec ~ 36.14t, and is well described by our strong-coupling 
analysis, which predicts an exciton with a size Chd ~ 2.78 
as the lowest optically excited state. However, one clearly 
sees in Fig. ||that Chdix) remains finite for large x. This 
wrongly suggests that the optically generated hole and 
double occupancy are not bound in this excited state. 
Similar problems arise with a density-density correlation 
function. Taking the difference between the correlation 
functions for an excited state and for the ground state, 
as in Ref. |l^, does not provide better results. 

A suitable quantity for our analysis is the correlation 
function for electron-hole excitations 



Cchix) = in Pi^i+, + {-lT\Pi+,,i 



(50) 



where |0) is the ground state, 
under investigation, and 



Pr 



is the excited state 



(51) 



creates an electron at site i and a hole at site j. Ob- 
viously, Coh {x) evaluates the importance of an electron- 
hole excitation with distance x in the excited state \n). 
This approach has already been used to study the struc- 
ture of excited states in semiempirical calculations for 
ladder-type PPP oligomersES. Here, we have calculated 
this correlation function using the symmetrized DMRG 
method to analyze the structure of excited states in the 
lattice model (^. 

An average electron-hole distance Cch can be defined 
by substituting Cch(a;) for Chd (a;) in Eq. (E^. We have 
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found that this method predicts correctly whether a hole- 
double occupancy pair is bound or unbound in the limit 
of large Mott gaps {U » t,Vi,V2). The advantage of 
the correlation function Cch{x) over Chd{x) (or a density- 
density correlation function) becomes obvious in a system 
with a Gnite Mott gap. In Fig. || one sees that Cch{x) 
decreases exponentially and thus allows us to identify 
a bound excitation, while the correlation function be- 
tween hole and double occupancy Chd(a;) is dominated 
by the finite ground-state density of holes and double oc- 
cupancies. Therefore, we think that the correlation func- 
tion for electron- hole excitations Cch{x) and the anal- 
ysis of the average electron-hole distance Ceh provide a 
reliable approach to distinguish an exciton from an un- 
bound particle-hole excitation in correlated systems. In 
any case, this approach is more reliable than schemes 
based on the correlation function for hole and double oc- 
cupancy or a density-density correlation function. 

However, this approach cannot be used to determine 
the size of a bound pair accurately because Ceh(a^) is af- 
fected by very strong short-range fluctuations due to spin 
correlations and lattice and chain-end effects. [In our fig- 
ures we show a running average of Ceh{x).] Instead, we 
evaluate the exciton size Ccx as the average electron-hole 
distance [as in Eq. (47)] calculated for the exponential 
fimction which gives the best fit to Cch{x) 



Cch{x) cx exp(-A|a;|) , 

c - I 

1 — exp(— A) 



(52a) 
(52b) 



In the limit of large Mott gaps we have found that this 
approach reproduces the exact result (49b) for the exci- 
ton size, i.e., (^ex ~ Chd- Furthermore, applying this anal- 
ysis to systems with large but finite gaps yields results 
in agreement with our strong-coupling analysis. For in- 
stance, using the data for U = 40t and Vi = 2.5t in Fig. |^ 
our analysis gives Cex = 2.99. This value agrees within 
10% with the analytical result (49b) for U/t — > oo and 
Vi = 2.5t, Chd = 2.78. The difference can be explained as 
a correction of the order At/U. Therefore, we think that 
the analysis (|5^) of the correlation function for electron- 
hole excitations allows one to determine reliably the size 
of an exciton in a correlated system. 

We have calculated the average electron-hole distance 
Ceh of the lowest optically excited states for various values 
of the parameters U, Vi, and V2 in the regime of interme- 
diate Mott gaps. The analysis of ^ch yields predictions 
about the presence of bound particle-hole pairs which al- 
ways agree with the predictions based on the analysis 
of the excited-state energies. For an excited state above 
the absorption band threshold we have found that ^ch 
diverges with system size, confirming that this excitation 
is made of an unbound particle-hole pair. A representa- 
tive example for the various possible shapes of Cchix) is 
shown in Fig. ^j. For an unbound pair one clearly sees that 
Coh{x) remains finite for very large x of the order of the 
system size L. Cchix) only goes to zero as x approaches 



the system size (L = 256 in this example) because of 
chain-end effects. 

In contrast, when an excited state lies below the charge 
gap Ec, we have found that Cch remains finite for an infi- 
nite system, confirming that this excited state is made of 
a bound particle- hole pair (exciton). In this case, Cchix) 
decreases exponentially and we can use the analysis d52) 
to determine the exciton size Cex- A typical example is 
also shown in Fig. We note that both examples in 
Fig. H correspond to systems with a gap Ec « 0.66t, for 
which the field theory of Sec. IV is still valid and supports 



our identification of bound and unbound excitations. 

The exciton sizes Cox observed range from slightly more 
than one lattice spacing ag to about 20ao. As the binding 
energy SE, the exciton size Ccx does not depend in a sim- 
ple way on the Coulomb interaction strength and even 
displays a non-monotonic behavior as a function of Vi, 
first decreasing as Vi increases, then sharply increasing 
as Vi approaches the critical regime U — 2Vi + 2V2 ~ 0. 
As expected, Cex becomes larger and seems to diverge as 
5E vanishes, while it diminishes when SE becomes larger 
if the gap Ec is kept (approximately) constant. 

We note that the operator used in (|5^) is antisymmet- 
ric with respect to charge conjugation, as is the current 
operator (||). Thus, Ceh{x) is non-zero only for those ex- 
cited states \n) (including the excited states contributing 
to the linear optical conductivity) whose parity Pc un- 
der charge conjugation is opposite to that of the ground 
state. When a minus sign is substituted for the -I- sign 
in (^o|), Coh(a;) is non-zero for excited states which have 
the same parity under charge conjugation as the ground 
state. One can thus extend this scheme to excited states 
which contribute to the non-linear optical properties of a 
system. 




FIG. 6. Correlation function Cch{x) of the first optically 
excited state on a 256-site lattice for two typical cases: an 
exciton (solid) for U = 4t, Vi — 2.75t, and V2 = 1.5t, and an 
unbound particle-hole pair (dashed) for U = 3.5t, Vi = 1.4t, 
and V2 = 0. The center of the pair is in the middle of the 
lattice. 

In the field-theory limit, a measure of the size of an ex- 
citon can be obtained in the following way. For asymp- 
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totically large separation between particles, the holon- 
antiholon wave function has the form 

^xi^x2{xi,X2) = exp(ipia;i + ip2X2) 

+Sn{pi,P2) eyip{ip2Xi +ipiX2) , 

*2:i»2:2(2;i:a;2) = S't(pi,P2) exp(jpixi +^^22:2) , (53) 

where S'r^t are two-particle S'-matrix elements corre- 
sponding to reflection and transmission respectively. In 
terms of rapidity variables defined by pj = M smh6j /vc, 
the S'-matrix only depends on the difference 0i2 — — ^2 
and has poles at 



i7r(l — n^) 



(54) 



An exciton state with energy Af„ cosh 9 and momentum 
M„sinh6'/?;c is obtained by choosing 



pi = — smh(t/ -|- I ^ 



P2 = — smh^fc* - I , 

Vc 2 



(55) 



Inserting these values in the wave function yields an ex- 
ponential decay in \x2 — xi\ with correlation length 

Vc 



FT 



Mcos 7r</2 
2vc 

V^' - Ml ■ 



(56a) 
(56b) 



We see that this size diverges when we approach from be- 
low the coupling constant ^, at which the n-th exciton is 
first formed. For example, the first exciton splits off the 
holon-antiholon continuum at ^ = 1 and its size displays 
a square-root divergence for ^ ^ 1 



croc 



1 



Ml ^ A" 



(57) 



A -Ml 

The non-monotonic behavior of the exciton size as a 
function of the Coulomb interaction strengt h is a lso ob- 
served in the small-gap regime. Equation (56a) shows 
that will decrease if the diminution of ^ dominates 
when the lattice parameter Vi increases, but becomes 
greater if the increase of the velocity Vc and the reduc- 
tion of the gap A prevail (see Sec. IV). 



To compare the field-theoretical predictions with our 
numerical method based on the correlation function 
Ceh(a;) we have numerically calculated the exciton en- 
ergy and size Ccx for several values of the parameters 
[/, Vi, and V2, corresponding to gaps A = £'c of the or- 
der of 0.6t. The exciton sizes range from lOao to 20ao. 
We have found that in all cases our numerical results ful- 



fill the field-theory relation ( 56b ) within 15% if we use 
Vc = 2tao as the charge velocity. This good agreement 
shows that our definitions ( |5^ ) and (56) are mutually 
consistent in the regime of small Mott gaps. Moreover, 
the numerical method yields reasonable estimates of the 
exciton size even in the regime where such calculations 
become laborious. 



VI. CONCLUSIONS 



We have investigated excitons in the optical conduc- 
tivity spectrum of one-dimensional Mott insulators us- 
ing two new reliable methods, the dynamical density- 
matrix renormalization group and the field-theoretical 
form-factor bootstrap approach, supplemented by two es- 
tablished techniques, a strong-coupling analysis and the 
symmetrized DMRG. Mott-Hubbard excitons can be un- 
derstood with the simple picture of a bound pair of spin- 
less bosonic excitations with opposite charge, in analogy 
to the bound electron-hole pair of Wannier exciton the- 
ory. However, the properties of Mott-Hubbard excitons 
are not as simple as those of Wannier excitons because 
both the Mott gap and the force between charged ex- 
citations originate from the Coulomb repulsion between 
electrons and are thus interdependent. In particular, the 
smallness of an exciton binding energy is no indication 
for the strength of the Coulomb repulsion in a Mott in- 
sulator. 

We may compare our results to experiments for poly- 
diacetylene chains in their monomer matrix. For an es- 
timate of the exciton size we use the field-theoretical re- 
sult (^) which is fairly independent of the details of the 
lattice structure and the range of the electron-electron 
interaction. We put Vc ~ ^ff— P 2tao with a 7r-electron 
band width = 4t « lOeVEl, A = Wph = 2.4 eV, and 
Ml = Wex = 1.9 eV from experiment for 3BCMU-PDA0. 
The field-theoretical prediction for the exciton size. 



FT 



W 



Go 



7an « 9 A 



(58) 



is in fair agreement with the experimental value 12 A if 
we put ao = 1.3 A for the average bond length between 
the carbon atoms on the PDA chain. The difference can 
be attributed to the limitations of this simple calculation 
or the uncertainty of the order of 20% for the experimen- 
tal valuescj. We therefore conclude that our many-body 
approach can be applied successfully to real polymers. 

Although the particular model studied here is too sim- 
ple to describe conjugated polymers accurately, we be- 
lieve that the concept of Mott-Hubbard excitons is rel- 
evant for these materials as the electron-electron inter- 
action significantly contributes to the formation of their 
optical gap. In any case, this approach is more realistic 
than the oversimplified Wannier theory and other sim- 
ple approaches that neglect or minimize the role of elec- 
tronic correlations in conjugated polymers. The many- 
body methods used in this work can be applied (and in 
part, have already been applied) to more realistic mod- 
els taking into account the polymer geometrical structure 
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and the electron-phonon interaction, and possibly addi- 
tional perturbations such as interchain couplings. We 
think that the optical properties of conjugated polymers 
will be successfully investigated using a combination of 
these many-body methods. 
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